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Abstract 

The nonlinear wave equation utt — c{u){c{u)ux)x = 0 determines a flow of conservative 
solutions taking values in the space i/^(]R). However, this flow is not continuous w.r.t. the 
natural distance. Aim of this paper is to construct a new metric which renders the 
flow uniformly Lipschitz continuous on bounded subsets of For this purpose, is 

given the structure of a Finsler manifold, where the norm of tangent vectors is defined in 
terms of an optimal transportation problem. For paths of piecewise smooth solutions, one 
can carefully estimate how the weighted length grows in time. By the generic regularity 
result proved in [7], these piecewise regular paths are dense and can be used to construct 
a geodesic distance with the desired Lipschitz property. 


1 Introduction 

Aim of this paper is to understand the continuous dependence of solutions to the nonlinear 
wave equation 

Utt - c{u)[c{u)ux)^ = 0. (1.1) 

Roughly speaking, the analysis in [8, 17, 22] shows that conservative solutions are unique, 
globally defined, and yield a flow on the space of couples {u,ut) G x L^(]R). For each 

conservative solution, the total energy 

F(t) = Jluf + c^(u)u^] dx (1.2) 

remains constant in time. Precise results in this direction will be recalled in Section 2. On the 
other hand, these solutions do not depend continuously on the initial data, w.r.t. the distance 
in the normed space 77^ x L^. 
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In the present paper we construct a new distance functional which renders Lipschitz continuous 
the flow generated by (1.1). We recall that, for solutions of the Hunter-Saxton or the Camassa- 
Holm equation, a similar task was achieved in [10, 13, 14, 20, 21]. 

Developing ideas in [13], our distance will be determined by the minimum cost to transport an 
energy measure from one solution to the other. While all previous papers dealt with first order 
equations, to define a suitable transportation distance between two solutions u,u of ( 1 . 1 ) one 
now faces three main difficulties: 

• At any given time t, each solution determines two distinct measures. These account 
for the energy of forward moving waves and the energy of backward moving 
waves. The distance between u{t) and u{t) should be measured by the minimum cost 
for transporting to /!(,_ and f/_ to /x(_. 

• The above double transportation problem is considerably complicated by the fact that, 
while the total energy is conserved, some energy can be transferred from forward to 
backward moving waves, or viceversa. These source terms must be accounted for, when 
designing an “optimal double transportation plan”. 

• As a wave front crosses waves of the opposite family, its speed can change. As a conse¬ 
quence, the distance between two corresponding fronts in u and u may quickly increase, 
making the optimal transportation plan more costly. To compensate for this effect, one 
needs to insert a weight function, accounting for the total energy of approaching waves. 

In Section 3 we introduce a Finsler norm on tangent vectors, related to an energy transporta¬ 
tion cost. Given a smooth path 7 : 0 i—)• (n®,nf), one can then define its weighted length 
jjyjj by integrating the norm of the tangent vector d'y/dO. Proposition 1, stated in Section 3 
and proved in Section 4, contains the key estimate, describing how the norm of a tangent 
vector grows in time. Assuming that, for 9 E [0, 1] and t G [0, T], all solutions u^(t, •) remain 
sufficiently regular so that the length of the path 7 * : 0 1 —)■ {u^{t), uf {t)) can still be computed, 
we obtain the bound 


hi < Ct 117°II, for all tG [0,T]. (1.3) 

Here the constant Ct depends only on T and on a bound on the x norm of the initial 
data. At this stage, it is natural to define the geodesic distance 

d*((u,ut), (u,ux)) = inf|ll 7 ll; j : [0,1] x , ^{0) = {u,ut), 7 ( 1 ) = (u, ut)|. 

(1.4) 

By (1.3) we thus expect that, for any two solutions of (1.1) and any t G [0,T], this distance 
should satisfy 

d*({u{t),ut{t)), {u{t),ut{t))J < Ct ■ d*(^{u{0),ut{0)) , (u(0), ut(0))^ . (1.5) 

This would imply that solutions depend Lipschitz continuously on the initial data, in the 
distance d*. 

To clinch this argument, one major difficulty must be overcome. Indeed, smooth solutions 
may well develop singularities in finite time, [19]. Given a path 7 ° of smooth initial data, 
there is no guarantee that at any time t G [0,T] the path 7 * will be regular enough so that the 
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tangent vectors d'y^/dO are meaningfully defined (see Fig. 1). We remark that a similar issue 
was encountered in the analysis of hyperbolic conservation laws [ 6 ]. For a path of piecewise 
smooth solutions with finitely many shocks, a weighted norm on a suitable family of tangent 
vectors was introduced in [5]. However, a lengthy effort was later required [9, 12], in order 
to construct paths of approximate solutions which retained enough regularity, so that their 
length could still be estimated in terms of these tangent vectors. 




Figure 1: Left: due to singularity formation, a smooth path of initial data 7 ° : 0 !->■ m®(0) may 
lose regularity at a later time T. In this case, the weighted length || 7 ^|| can no longer be computed 
by integrating the norm of a tangent vector. Right: by a small perturbation of the initial data, one 
obtains a path of solutions 0 1 —>■ which remain piecewise smooth, for all except finitely many values 
of 0 € [0,1] . 

In the present context, we can take advantage of the generic regularity results recently proved 
in [7]. These can be summarized as follows. 

(i) For an open dense set of initial data 

M^r)) = {uo,ui) e X (c^{R)nL‘^{R)'^ (1.6) 

the corresponding solution u = u{t, x) of (1.1) is piecewise smooth in the t-x plane, with 
singularities occurring along a finite set of smooth curves. 

(ii) Every path of initial data 9 1 —)• 7^(0) = {uq,u\) can be approximated by a second 
path 6 1-7 7 ^( 0 ) = ('UQj'uf) such that, for all but finitely many values of 0 € [0,1], the 
corresponding solution remains piecewise smooth on the domain [0, T] x R. 

Using this dense set of piecewise regular paths, we can thus define a geodesic distance on the 
space x L^, with the desired Lipschitz property. Our main results are contained in 

• Proposition 1, which establishes the basic estimate (3.22) on the size of tangent vectors. 

• Theorem 5, providing the bound (6.3) on how the length of a path of solutions can grow 
in time. 

• Theorem 7, showing that, by (7.6), the flow generated by the wave equation (1.1) is 
Lipschitz continuous w.r.t. the geodesic distance d*. 

We remark that, for hyperbolic conservation laws, the distance constructed in [5, 9, 12] is 
equivalent to the distance. On the contrary, our new metric is not equivalent to the norm 
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distance on H 1 xL2 . The completion of x w.r.t. the geodesic distance includes a family 
of measures. This should not come as a surprise, since it was already observed in [17, 22] that 
conservative solutions can occasionally be measure-valued. 

In Section 7 we compare the geodesic distance (1.4) with more familiar distances found in the 
literature. In one direction, we show that 

d*({uo,Ui), iuo,Ul)) < C ■ (\\uo - UoWjjI + \\uo - Uo\\wi-i + \\UI - Ul\\-L 2 + \\ui - UiWj^i) , 


for some constant C. On the other hand, let /i and fl be the positive measures having densities 
respectively 

uf + c^{u)u^ and uf + c‘^{u)u^ (1.7) 

w.r.t. Lebesgue measure. Then the geodesic distance d* dominates the Wasserstein distance 
between the two measures. Namely 


sup 



Cl ^ 


< 1 




( 1 . 8 ) 


All of the present analysis is concerned with conservative solutions to (1.1). For dissipative 
solutions, studied in [15, 19, 26, 27], the continuous dependence for general initial data in 
X remains an open question. For scalar conservation laws, an entirely different approach 
to continuous dependence, relying on an formulation, was developed in [2, 3, 4]. 


2 Conservative solutions to the nonlinear wave equation 


In this section we review the main results in [7, 8, 17] on the Cauchy problem for the quasilinear 
second order wave equation 

utt-c{u){c{u)ux)^ = 0, ( 2 . 1 ) 

with initial data 


u( 0 ,x) = uo(x), ut{ 0 ,x) = ui{x). ( 2 . 2 ) 

Here c : M i—)• M+ is a smooth, uniformly positive function, such that 


Consider the variables 


c{u) 

IV 

o 

o 

V 

o 

R = 

Ut + c{u)u 

S = 

Ut — c{u)u 


X ) 
X ) 


so that 


ut 


R + S 
2 


By (2.1), the variables R,S satisfy 


'^x 


R-S 

2 c 


Rt-cR. = 

St + cS, = ^^(S^-R^). 


(2.3) 

(2.4) 

(2.5) 


( 2 . 6 ) 
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Multiplying the first equation in (2.6) by R and the second one by S, one obtains balance laws 
for and 5^, namely 


[R^)t-{cR% = f^{R^S-RS^), 
{S% + {cS% = -f^{R^S-RS^). 


As a consequence, for smooth solutions the following quantity is conserved: 


2 I 2 2 

e = u-j. + c 


i?2 + 52 
2 


(2.7) 


( 2 . 8 ) 


We think of i?^/2 and 5^/2 as the energy of backward and forward moving waves, respectively. 
These are not separately conserved. Indeed, by (2.7) energy is transferred from forward to 
backward waves, and viceversa. The main results on the existence of solutions to the Cauchy 
problem can be summarized as follows. 


Theorem 1. Let c : M M be a smooth function satisfying (2.3). Assume that the initial 
data uq in (2.2) is absolutely continuous, and that {uq)x £ , ui G L^. Then the Cauchy 

problem (2.1)-(2.2) admits a weak solution u = u{t,x), defined for all (t,x) G MxM. In the t-x 
plane, the function u is locally Holder continuous with exponent 1/2. This solution 1u{t, •) 
is continuously differentiable as a map with values in for all 1 < p < 2. Moreover, it is 
Lipschitz continuous w.r.t. the distance, i.e. 

||u(f, •)-m(s,-)||l 2 < L\t-s\ (2.9) 

for all t, s G M. The equation (2.1) is satisfied in distributional sense, i.e. 



(c(u)(/>)^c(u) Ux 


dxdt = 0 


( 2 . 10 ) 


for all test functions cj) £ C). The maps t ut{t, •) and 1Ux{t, •) are continuous with values 
in L/^^(M), for every p £ [1, 2[. 


Theorem 2. In the same setting as Theorem 1, a unique solution u = u{t,x) exists which is 
conservative in the following sense. 

There exists two families of positive Radon measures on the real line: and de¬ 

pending continuously on t in the weak topology of measures, with the following properties. 

(i) At every time t one has 

/ OO 

-oo 

(a) For each t, the absolutely continuous parts of and //(,_ w.r.t. the Lebesgue measure 
have densities respectively given by 

R^ = [ut + c{u)uxf‘, 5 ^ = {ut - c{u)uxf‘. ( 2 . 12 ) 

(Hi) For almost every t G M, the singular parts of /xL and are concentrated on the set 
where c'(tt) = 0. 


ulix) + {c{uo{x))uo,xix))‘ 


dx. 


( 2 . 11 ) 
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(iv) The measures fit. and /i^,_ provide measure-valued solutions respectively to the balance 
laws 


= f,{R^S-RS^), 

'nt + icr])x = -^{R‘^S - RS^). 


(2.13) 


The existence part of the above theorems was proved in [17]. The uniqueness of conservative 
solutions has been recently established in [8]. 


Remark 1. By (2.13) the total energy, represented by the positive measure /r* = /i(|_ + is 
conserved in time. Occasionally, some of this energy is concentrated on a set of measure zero. 
At a time r when this happens, has a non-trivial singular part and hence its absolutely 
continuous part satisfies 


j u‘fir,x) + c^{u{T,x))uliT,x) 


dx < Eq . 


The condition (iii) puts some restrictions on the set of such times r. In particular, if c'{u) ^ 0 
for all u, then this set has measure zero. 


Remark 2. For any t > 0, the conservation of the total energy implies 

\\Mt)\\h < Eq = j{ul + c^{uo)ul ,^)dx. (2.14) 

Hence (2.9) holds with Lipschitz constant L = y/Eo- Moreover, one has the bounds 

\\uit,-)\\-L2 < \\uo\\i^2 + t y/^, \\Ux{t, ■)\\-L^ < . (2.15) 

Co 

This yields an a priori bound on ||u(t, •)||j 7 i, and hence on ||u(t, •)||loo, depending only on time 
and on the total energy Eq. In turn, since the wave speed c(-) is smooth, we obtain an a priori 
bound on c{u) and |c'(tt)|. 


3 First order variations 


For simplicity, in this section we consider solutions of (2.1) with bounded support. More 
precisely, we shall assume that all our solutions satisfy 

u{t,x) = 0 for a; ^ [0 ,Lq], t€[0,r]. (3.1) 


Because of finite propagation speed, this is hardly a restriction. 

Let {u,R,S) provide a smooth solution to (2.1), (2.4), and consider a family of perturbed 
solutions of the form 


= u-\- sv o(e ), 


f R^ = R-h sr -h o(e), 
\ S'^ = 5" + es + o(e). 


(3.2) 
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From (2.5) it follows 


+ R + S r + s 


^ -S^ R-S r-s R-S ,, , 

2 c(ri^) 2c{u) 2c(u) 2c2(u) ^ ^ ^ ^ 

Under the assumption (3.1), given r,s, the perturbation v is uniquely determined by 

(R — S)c'(u) r — s 


Vt. = - 


Furthermore, we have 


2c^{u) ^ 2c(ri) ’ 

r + s 


v{t,0) = 0. 


(3.3) 

(3.4) 

(3.5) 

(3.6) 


A direct computation shows that the first order perturbations u, s, r satisfy the linear equations 


C Vxx — 2cC Ux'^x H“ 2cC Uxx^^ • 

n - c{u)rx = c'RxV + - S'^)v + |^(^^ - 5's), 


St + c{u)sx = - c'SxV + 




(3.7) 


(3.8) 


By the assumptions (2.3) on the wave speed c{u), all functions c'/4c, c"/4c, (c')^/4c^, are 
smooth functions of u. 


We shall introduce a weighted norm on tangent vectors r, s, which takes into account the total 
energy of waves which are approaching a given wave located at x. This is described by the 
weights 


r* + C50 


W (x) = 1 + / S‘^{y)dy, >V^(x) = 1+1 R^{y)dy. 

J —OO J X 

In addition, consider the function 

\c'\\R^S - S‘^R\ , , 

a{t) = / --- -{t,x) dx. 

J-oo 2 c 

As proved in [8], the function 


r I—)• 


rr r+co 

'0 J-oo 


2c 


{R^S - RS^ 


{t, x) dx dt 


(3.9) 


(3.10) 


is Holder continuous and absolutely continuous on bounded time intervals, and has sub-linear 
growth. In particular (see (3.11)-(3.12) in the proof of Lemma 1 in [8]), one has 


ra{t) 

Jo 


dt < Ct , 


(3.11) 
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for some constant Ct depending only on T and on the total energy Eq. By (2.7) it follows 


wr - cw- 
Wt + cW+ 


- 2cS^ + 

- 2cR^ + 


r ^{S^R-R^S)dy < 

J —oo 

^+oo / 

J -{R^S-S^R)dy < 


— 2coiS^ + ait), 


— 2cqR? + a{t). 


(3.12) 


On the space of tangent vectors {v, r, s) we introduce a Finsler norm, having the form 


(u, r, s) 


= _ inf 

{u,R,S) r,s,w,z 


{r,w,s,z) 


iu,R,S) 


(3.13) 


where the infimum is taken over the set of vertical displacements f, s and shifts w, z which 
satisfy 

r = r - wRx +- z)S'^ , 

(3.14) 

S = S — zSx +-^{w — z)R‘^ . 

This norm is defined as 


{r,w,s,z) 


{u,R,S) 

Ki h\w\{l + R‘^)W-+ \z\{l + S^)W+^dx 


+K2 I {|f|W- + |s|W+} 


dx 


+R-3 


/ 


V + 


Rw Sz 


I (1 + R^) W + (1 + ^2) W+|dx 


+ K4 




W + 




W+ \dx 


+K5 j I Rwx + “ z)SR 


W + 


Szx + -^{w - z)RS 
4c^ 


W+ Ux 


+ K6 


2Rr + R?Wx + -^R?S{w — z) 
4c^ 


>V" + 


2Ss + S^zx + ^S^Riw - z] 
4c^ 


>V+ \dx 


— Rih + R2I2 + R3I3 + ) 

(3.15) 

for suitable constants ki, ..., Ke to be determined later. 

To motivate (3.13), consider a profile R and a perturbation R^, as shown in figure 2. In first 
approximation, R^ ^ R + er. Notice that we could also obtain the profile R^ starting from 
the graph of R, performing a horizontal shift in the amount ew and then a vertical shift in 
the amount sf, provided that 

r = f — wRx ■ (3.16) 
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As a first guess, one could thus define a norm ||r||-j- by optimizing the choice of r,w, subject 
to (3.16). However, a detailed analysis has shown that this approach does not work. Indeed, 
it does not take into account the fact that, when backward and forward moving waves cross 
each other, by (2.6) their sizes R,S are modihed. Compared with (3.16), the additional term 
in the hrst equation of (3.14) accounts for this interaction. Notice that tc — z is the relative 
shift of backward w.r.t. forward waves. 



Figure 2: A perturbation of the i?-component of the solution to the variational wave equation. 
We now explain the meaning of each integral on the right hand side of (3.15). 

• The integral of |tc|(l + A^) can be interpreted as the cost for transporting the base 
measure with density 1 + E? from the point x to the point x + ew{x). 

Similarly, the integral of \z\{l + A^) accounts for the cost of transporting the measure 
with density 1 + 5^ from x to x + ez{x). 

Here, as in all other terms, we insert the weights coming from the interaction 
potential. 


• I 2 accounts for the vertical shifts in the graphs of R, S. We interpret the integrand as 
the change in arctan A times the density (1 + R^) of the base measure. Notice that here 
the factor (1 + i?^) cancels out with the derivative of the arctangent. 


I 3 accounts for the changes in u. Observe that 


e ^[u‘^{x + ew{x)) — u{x)] « v{x) + Ux{x)w{x) = v{x) + 


R{x) — S{x) 
2 c(u(x)) 


w(x). 


This can be written in the form 


V + 


R-S 


2 c 


R 


^ = "• + 2c'“ - 2c^ ' + 


S \ S{z — w) 


2c 


(3.17) 


Notice that the last term on the right hand side of (3.17) does not appear in I 3 . In 
fact, the last term relative shift term coming from the equation (2.4). 

Subsequent computations will show that this term is inessential, because its contribution 
can be bounded by the decrease in the interaction potential. In an entirely similar way 
we obtain 


e {x + ez{x)) — u{x)\ « v{x) + Ux{x)z{x) 


v{x) + 


R{x) — S{x) 
2 c{u{x)) 


z{x), 
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R-S 


R S 

V H- w -z 

2 c 2c 


+ 


R{z — w) 
2 ^ 


• Iq accounts for the change in base measure with densities R^ and S'^, produced by the 
shifts 'w,z. To see this, assume that the mass with density R^ is transported from x to 
X + ew{x). If the mass were conserved, the new density should be 

{R‘^)‘^ {x + £w{x)) = R^{x) — £w{x){R?)x{x) — £Wx{x)R?{x) + o{£). (3.18) 


In addition, if the mass with density is transported from x to x + £z{x), by (2.7) the 
crossing between forward and backward waves yields the source term 

^^(R^S-RS^)-^. (3.19) 

On the other hand, if we shift the graph of R horizontally by £w and then vertically by 
ef, the new density will be 

{R'^)‘^{x + £w{x)) = R^{x) — £w{x){R‘^)x{x) + 2£R{x)r{x) + o{£). (3.20) 


Subtracting (3.18)-(3.19) from (3.20) we obtain the expression 



(3.21) 


• The integrals R and R does not seem to have a clear geometric interpretation. R is 
somewhat related to the change in Lebesgue measure produced by the shifts w, z, while 
R is related to the change in base measure with densities R and S, produced by the 
shifts w, z. As shown by our subsequent computations, these two additional terms must 
be included in the definition (3.15), in order to estimate the time derivatives of R and 

R- 


Our goal is to prove 

Proposition 1. Let {u,R,S) be a smooth solution to (2.1) and (2.6), and assume that the 
first order perturbations {v,r,s) satisfy the corresponding linear equations (3.7)-(3.8). Then 
for any r > 0 one has 


(u(T),r(T),s(r)) 


< exp<^C'r+ / a(s)ds|>- (u(0), r(0), s(0)) 

(u(t),R(t),S(t)) I Jq 


with a constant C depending only on the total energy. 


{u{o),Rio),sm 

(3.22) 


Toward the proof, the main argument goes as follows. At time t = 0 let a tangent vector 
(u(0), r(0), s(0)) be given. By the definition (3.13), for any e > 0 we can find shifts wq, zq and 
perturbations fo, sq satisfying 


(fo,u;o,so,^o) 


(«(0),i?(0),5(0)) 


< e + 


(u(0), r(0), s(0)) 


(«(0),K(0),S(0)) 


(3.23) 
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together with the constraints 

f ^(0) 


= fo - woRx{0) + ^{wo - zo)5'^(0), 


(3.24) 


[ s(0) = So - ^;o*S'a;(0) + - ^;o).R^(0). 

In order to prove (3.22), for any t G [0, r] it suffices to hnd shifts w^t), z{t), together with f{t), 
s{t) satisfying (3.14) and the initial condition (3.24), so that 


A 

dt 




< {C + a{t))- {r{t),w(t),s{t),z{t)) 


(3.25) 

These shifts w(t),z{t) will be obtained by propagating along characteristics the shifts wo,zo 
in the initial data. More precisely, we choose w, z to be the solutions of the linearized system 


wt - c{u)wx = - c'{u) {v + UxW) , 
Zt + c{u)Zx = c'{u){v + UxZ) , 


(3.26) 


with initial data 


w{0, x) 


wo{x ), 


(3.27) 


[ z{{),x) = zo{x). 

By (3.8) and the identities (3.14), this determines the evolution equation for r,s. 

In the next section, by carefully estimating the time derivatives of all terms in (3.15), we shall 
prove that (3.25) holds. In turn, this will yield (3.22). 


4 Estimates on the norm of tangent vectors 


The first part of the proof of (3.25) is largely computational. Using the evolution equations 
(2.1), (2.4), (2.6) for u,R,S, and (3.8), (3.26) for r,s,w,z, together with the identities (3.14), 
we estimate the time derivative of each integral in (3.15). 

1. To estimate the time derivative of R (shift in the base measure), using (3.26) we first 
compute 

(u;(l + i?2))^ _ (cu;(l + 

= (wt - CWx){l + i?^) + w [{R‘^)t - {cR‘^)x] - WCx 

= -c'(v + ?-^w) (1 + i?2) + —w{R^S -RS'^ -R + S) 

V 2c / 2c 

= -d{v + —w- —. 2 ) (1 + R^) + —w{2R^S -RS"^ -R + 2S) - —zS{l + R^). 

\ 2c 2c J 2c 2c 
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Thanks to (3.12) we obtain 


dt 


I |u;| (1 + R^)W- dx < 0(1) • J |w;| (l + \R^S\ + + \R\ + |S|)W- dx 

+0(1) • J |z| {\S\ + \R^S\)W+ dx + 0(1) • J 


Rw Sz 

V ^- 

2c 2c 


(1 + R^)W- dx 


+a{t) J |t(;| (l + W dx — 2co J |i(;| (l + dx . 


(4.1) 


2. To estimate the time derivative of I 2 (change in arctan), using (3.8) we first compute 
(r + wR^)^ - (c(r + wR^^))^ 

= [rt - (cr)^;] + {wt - cWx)Rx + w[{Rx)t - {cRx)x] 


-c' u + 


R-S 


2c 


w ] Rx + w 


c"c - (c')^ R-S 


(i?2 _ 52) 


4c 


{2RRx - 2SSx 




/ E? _ Q 

^r2_s^)(v + ^^w 


2c 


Next, 


(4.2) 


8 c 2 


{w-z)S‘^ ] - ( c^{w-z)S 


c"c - 2(c') 


t 

.'\2 


- {ut - CUx){w - z)S‘^ + -^{wt - CWx)S‘^ “ 


+ (cS'^),^] - ^2{w - z){cS‘^)^ 


8c2 


8c2 


/\2 




8 c^ 


8c2 


V + 


R-S 


2c 


w\S^- 


r,l\2 


8 c 2 


R-S 


2c 


z S^ 


+^/xS^ - - ^){R^S - RS^) - ^(w - z)(RS^ - S^} - ^(w - z)SSx 


/^2 


(cO 


/^2 


(4.3) 
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By (3.14), combining (4.2) with (4.3) we obtain 


h - {cr)x 


= [{r+ wR^)^-{c{r+ wRx))^] - 




c „ . c , c"c— (c ')^.„2 c 2 ' ^ R — S 


= -Sv,{R, - SJ + -SXr -,)+ 


-(R^-S^) U + 


2c 


w 


c"c-2{c'f, ^-'^2 

^ -{w- z)S-^ + 


8c3 


8 c 2 


2c 


3 (c')Vn R-S. ^2 c' ^2 

' ‘ 2v -\ - (w + z) ) 52 — —ZxS^ 

/ 4 ^ 


(c') 


/^2 


+^(?/; - ^)(i?'5 - RS^) + ^(w - z)(RS^ - S^) + g(w - z)SSx 
= ^{wSRx - zSSx) - + Y^S ({r - s) - (wR^ - zS^) + - z){S‘^ - R^) 


,'^2 


Sc^ 




{w - z)S^ 


(c')2 / R-S, , > 2 




w — z, 


c"c - (c-)^ , / jw _ & \ _ c"c - 2(cy g 2 / ^ _ S£\ 

4c2 V 2c 2c J 4c2 V 2c 2c J 

+0{l)-{\w\ + \z\){l + \R^S\ + \RS^\). 


We thus conclude 


dt 


j\r\W-dx = 0{l)- j\Sr\W-+ 0(1) ■ j 


2Ss + S^Zx + —S^R(w-z) 
4c 


(4.4) 

W+dx 


+0(1) • I 


52 


Rw Sz 

V H- 

2c 2c 


W+ dx + 0{l)- J 


R^ 


Rw Sz 

V H- 

2c 2c 


>V"dx 


+0(1) • j\w \(1 + |i?25| + \RS^\) W dx + 0(1) • j |z| (1 + \R^S\ + \RS^\) W+ dx 

+a{t) J |f|>V“(ix — 2co J |f| S^W~dx. 


13 
























3. To estimate the time derivative of (change in u), using the identities in (3.5)-(3.6) for vt 
and Vx, we hrst compute 


vt-cvx = s + —{R- S)v. 


Next, by (2.4) and (3.26) we obtain 


(Rw Sz 


Rw Sz\ 


(4.6) 


= ^w{Rt - cRx) - ^z{St + cSx) + zSx 
2c 2c 


R S c' 

+ - CWx) - —{zt + CZx) + Szx - -^{RSw - s‘^z) 


(4.7) 


= -^w{R^ - S^) - ^z{S^ - R^) + zSx 


c'( R-S \ 
7?(vH-;-rcl - 


) ■ yA "+ 


2 c 


Finally, by (2.6) it follows 


(1 + R^t - ic(l + R^h = ^(R^S-RSA-^(R-S). 

2c 2c 

Putting together (4.6)-(4.8) and using (3.14) one obtains 


V + 


Rw Sz 
2c 2c 


){i + rA 

' J t 


civ + 


Rw Sz 
2c 2c 


+ R^ 


Vt - CVx + 

Rw Sz 


Rw Sz 
2c 2c 


l-i 


Rw Sz 
2c 2c 


(i + rA 


+ u + 


2 c 2c 


{1 + R%-{c{1 + RA), 


s + ^v{R -S) + ^w{R^ - SA - - RA + zSx 

2c 8c^ Sc‘^ 


-—R{ V + 


R-S 




2 c 


+ uT 


2 c 

Rw Sz 
2c 2c 


-w — 


{1 + RA 


^^{r^S-RSA-y^{R-S) 


i+^{z- m)S^ -^-s{v + ^-^) + [Sz. + -^{w- z)RS) 


(1 + rA 


Rw Sz 


^^(R^S-RSA-Y^iR-S) 


(4.8) 
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We thus conclude 


dt 


V + 


Rw Sz 


(1 + R^) W dx 


< 


Rw Sz 


j |5|(1 + R^)W+ dx + 0(1) • J 

+0(1) ■ j\w\S^{l +R^)W-dx + 0{l) ■ j\z\S^{l +R^)W+ dx 


(1 + |/?| + |5| + |i? 25 | + \RS‘^\) n;- dx 


+ 0 ( 1 ) 

+a{t) [ 


Szx + - z)RS 

4c^ 


(1 + R^)W~ dx 


V + 


Rw Sz 
2c 2c 


{l + R^)W- -2co 


V + 


Rw Sz 
2c 2c 


(1 + R^)S^ W-dx . 


(4.9) 


4. To estimate the time derivative of 14 ,^ recalling (3.26) we first compute 




R-S 


2c 


w — c 


{R - S)c' r - s c' 2 

- T-^v H- -r^{R - SYw 

2c2 2c 4c3 ^ ^ 


(4.10) 


-^(^xW^ - SxW) - • 

Moreover, by (2.4) and (3.26), one has 


1 - 




,/'i 2 


(£l 

4c^ 


R-S 


4c^ 


-^zS - c-^zS 


4c^ 


(c') 


/^2 


M 2 


(cO 


V H - r:—w ] S - j^w{R - s ) - - S)wS - —wSx , 


2c "J ~ 16c3 

4c2i ^ 4 c2 


2 c 


(4.11) 




(c^i d d 

^-^(R - S)zS - zS^ - Z:,S . 

8c 3 ^ ^ 2c 2c 


(4.12) 
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Combining the identities (4.10)-(4.12) and recalling (3.14), we obtain 




- t 


c' _ c' . c' 


= —S — —f+—S(Wx + -^{w — Z 
2 c 2c 2c V 4c2^ 

j , J 


„/ Rw Sz 

i? cH- 

V 2c 2c 


(4.13) 


c"c - (c')2 


+ 


2 c2 

c"c - (c')2 


Rw Sz\ c"c —2(c')^^/ Rw Sz 


J \2 


4^3 - ^) - 


By the previous analysis, thanks to the uniform bounds (3.12) on the weights, we conclude 


dtj 




W- dx 


< 0(1) ■ j \f\W-dx + Oil)- j \s\W^dx + Oil) ■ j S w^ + ^{w - z)S 
+ 0 ( 1 ) • J Szx + “ z)RS W~^dx + 0 ( 1 ) • J Rwx + “ z)RS 

+o{iy I 


W-dx 


W-dx 


Rw Sz 

V H- 

2c 2c 


\R\W-dx + 0{iy 


Rw Sz 

V 4- 

2c 2c 


\S\W+dx 


+0(1) • j |u;| {\RS\ + 5^) W dx + 0(1) • j |^| (|i?5| + 5^) W+ dx 


+a{t) 


Wx + 


W dx — 2co • 


Wx + 


52 w dx. 


(4.14) 
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5. To estimate the time derivative of Is, using (4.13) we compute 
r(^Wx + - z)s'^ + Rc(wx + ^{w-z) 




+^r(^Sz^ + ^(w - z)SR'j - ^r(^Rw^ + ^(w - z)SR'j 

c"c — f c"c — 2(c')^ Rw Sz 

2c2 V 2c 2c) 2c2 V 2c 2c 


4c'^ Sc'* 




= -«»-:r 

2c 4c 


{2Rr + 


+^5(i?rc^ + -^{w - z)RS) + Yc^(Sz^ + -^{w - z)SR^ 

2c2 V 2c 2c) 2c2 V 2c 2c / 


(4.15) 


c"c - (c')^ „2 ^ ('^0^ e2 E>^ ^ c 2 / , C' , 

+—4^3—^ - ^) - r" + ~ ^ 


/ 

4c 


4c2 
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We thus conclude 


A 

dt 


Rwx + - z)RS 


W- dx 


< 0(1) • j |i2s| W dx + 0(1) ■ j 2Rr + R^w^ + -^{w - z)SR^ 
+0(1) • y Rwx +-^{w — z)RS \S\W~dx 
+0(1) • y Szx + -^{w — z)SR \R\W~dx 


W~dx 


+o{i)- 1 


Rw Sz 


R^W-dx + 


om.f 


Rw Sz 

V -\ - 

2c 2c 


\RS\W+dx 


+ 0 ( 1 ) • y (|u;| + \z\) (1 + i? 2 )(l + ^ 2 ) ^ . y ^ 

+a{t) ■ y 


Wdx 


s 

f 

c' 

Rwx + -r-^(w — zjRS 
4c^ 

W dx — 2cq ■ 

Rwx + -r-^{w — z)RS 
4c^ 


W dx 


(4.16) 
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6. Finally, to estimate the time derivative of Iq (change in base measure with density R^), we 
compute 


^2i?f + R^Wx + “ z)SR^ + + R?Wx + “ z)SR^^'^ 

= {Rt - cRx)(2r + Rwx + ^{w - z)SR^ 


+R 


2(j^t - (cr)x) + (^{Rwx)t - {cRwx)x'j + “ z)SR'j 

= ^^{R^-S^)(2r + Rwx + ^{w-z)SR) 


a 

t V4c 


{w — z)SR)j 


+- 


!/ /2 
C C — C 


2 c 2 


Rw Sz\ c''c — 2 c!^ 


R^v + -- 

2 c 2 c 


2 c 2 


2 c 2 c 


jr „ „/2 


^C-C-_ y;) _|_ 1- l-R^S‘^{w — z) 

4c-^ 8 c-^ 

-—r(2Ss + S^Zx + -^S^Riw - z)) + -SRr 
2c ' 4c ' c 


c' „ 2 ;- c' 


+ - ^)SR^) + y^Sr{Rwx + ^(w - z)RSJ 


H —R s — —R{ 2Rf + R 
2c 4c V 

+^^R^{Szx + ^iw-z)SR) 


c''c—{c'Y o/ Rw Sz\ d'c — 2 {d'f‘ Rw Sz 

^ -R^[v + - -— I +- ^ ^ SR^{v + - -— 


2 c 2 


2 c 2 c 


2 c 2 


2 c 2 c 


+ ^ ^ R^S{w - z) - ^^^S‘^R^{w - z) - ^S‘^r(wx + -^{w - z 

Act sc^ 4c V 4c^ 


+ - 


c"c - 2c'2 


2 c 2 


V + 


2 c 2 c ^ 2 c 2 c 


2 c 2 c 

j 


c2~ c' 

-5^r H-. 

2 c 2 c 


{2Rr + R^wx + -^R^Siw - z)) - ^R f2Ss + S^Zx + ^RS^{ 


2 c V 


2 c V 


4c2 


[w - z, 


(4.17) 
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This yields the estimate 


d_ 

dt 



z)SE? 


W~dx 


< 0{1)- 




Szx + 


W^dx 


+ 0 ( 1 ) 

+ 0 ( 1 ) 

+ 0 ( 1 ) 

+ 0 ( 1 ) 


W-dx 


Rw Sz 

V H- 

2c 2c 


J Rwx + — z)RS 

\R^S - RS^\W- dx 

J S^\f\ W dx + 0(1) • j R^\s\ W+ dx 
[ |5| 


2Rr + R^Wx + -—-^R^S(w — z) 
4c^ 


W-dx 


(4.18) 


+ 0 ( 1 ) 



255 + ^ 2 ^^^ 0 ^ 2 / 

4c2 




W+dx 



Figure 3; A graphical summary of all the a priori estimates. If a lower box Ik is connected to an 
upper box Ig, this means that the integral R is used in order to control the time derivative R = ^R- 
If £ € then ik and It are connected by a solid line. If £ € then Ik and It are connected by a 
dashed line. 
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7. We keep track of all the above estimates by the diagram in Fig. 3. Recalling (3.15), the 
weighted norm of a tangent vector can be written as 


{f,w,s,z) 


{u,R,S) 


— Rlh + ^2-^2 + R3I3 + RaIa + K5/5 + KqIq 


^ (/ Jk >V" dx + J J+ W+ dx^ 


(4.19) 


where are the various integrands. According to the estimates (4.1), (4.5), (4.9), (4.14), 

(4.16), and (4.18), the time derivative of each can be estimated as 

ik < + 1*^1^ W-dx + j J+ (1 + \R\) W+ dx'^ 




(4.20) 


+ ^ 0(1) + R^) W-dx + j J+ (1 + S^) W+ dx'] 

+a{t)Ik - 2 co (^J S^J^: W-dx + j R^J+ W+ dx^ . 

Here R14 C {1,2, ...,6 } are suitable sets of indices, illnstrated in Fig. 3. By direct in¬ 
spection, we see that the set-valued map k has no cycles. Indeed, the composition 

rIo R^o R^ yields the empty set. 

By choosing a constant <5 > 0 small enough, we thus obtain a weighted norm 


{r,w,s,z) 


{u,R,S) 


= h + 6I2 + S^l3 + dh + 6^ h + h 


(4.21) 


which satisfies the desired inequality (3.25). This completes the proof of Proposition 1. □ 


5 Tangent vectors in transformed coordinates 

Given any path 6 1 —?> u®, 0 G [0,1] of smooth solutions to (1.1), the analysis in the previous 
section has provided an estimate on how its weighted length increases in time. However, even 
for smooth initial data, it is well known that the quantities ut , Ux can blow up in finite time 
[19]. When this happens, a tangent vector may no longer exist; even if it does exist, it is not 
obvious that our earlier estimates should remain valid. Aim of this section is to address these 
issues. Roughly speaking, we claim that 

(i) Every path of solutions 9 ^ can be uniformly approximated by a second path 9 ^ 
such that, for all but finitely many values of 0 G [0,1], the solution is piecewise 
smooth, with “generic” singnlarities. 

(ii) If all solutions are piecewise smooth, with “generic” singularities along finitely many 
points and finitely many curves in the t-x plane, then the tangent vectors are still well 
defined and their norms can be estimated as before. 
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A precise formulation of (i) was recently proved by the authors in [7]. The proof is based on 
the representation of solutions to ( 1 . 1 ) in terms of a semilinear system with smooth coefficients 
[17], followed by an application of Thom’s transversality theorem. We review here this basic 
construction, and the characterization of generic (structurally stable) singularities [16]. 

To deal with possibly unbounded values of R,S in (2.4), following [17] it is convenient to 
introduce a new set of dependent variables: 

a = 2arctani?, (5 = 2arctanS'. (5.1) 

Using (2.6), we obtain the equations 


2 , , 

c' R^ - 

(5.2) 

c (yx — 1 -|- ^ ) — 

2 c l + i ?2 ’ 

2 . 

C' 52 - i?2 

(5.3) 

IJt + cjJx — ^ ~ 

1 + 

2 c 1 + 52 ■ 


We now perform a further change of independent variables. Consider the equations for the 
backward and forward characteristics: 


X = — c{u ), = c{u) , (5.4) 

where the upper dot denotes a derivative w.r.t. time. The characteristics passing through the 
point (t, x) will be denoted by 

s !-)■ x~(s,t,x}, s x~^{s,t,x), 

respectively. We shall use a set of coordinates (AT, Y ) on the t-x plane such that X is constant 
along backward characteristics and Y is constant along forward characteristics, namely 

( Xt- c{u)Xj, = 0 , 


(5.5) 


[ Yt + c(u)Yx = 0 . 

For example, one can dehne X, Y to be the intersections with the x-axis, of the characteristics 
through the point {t,x), i.e. 

X{t,x) = x~{0,t,x), Y{t,x) = —x~^{0,t,x). (5.6) 

More generally, one can consider strictly increasing functions x X(x) and x Y (x) and 

define 


X{t,x) = X(x (0,t,x)), Y{t,x) = y(—x'''(0,t,x)) . (5.7) 


For any smooth function /, using (5.5) one finds 

r ft + cU = fxXt + fYYt + cfxX^ + cfyY^ = {Xt + cX^)fx 

\ ft-cU = fxXt + fyYt-cfxX^-cfyY^ = {Yt-cY^)fy 
We now introduce the further variables 


V 


l + R^ 
X^ 


q 


1 + 52 

-Ux 


‘^cXxfx , 

(5.8) 

-2cY^fy . 


(5.9) 
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Notice that the above dehnitions imply 
p (1 + cos a)p 


1 


l + R^ 


-1 

y:. 


(1 + cos I3)q 


1 + 52 


(5.10) 


Starting with the nonlinear equation (2.1), using X,Y as independent variables one obtains a 
semilinear hyperbolic system with smooth coefficients for the variables u,a, I3,p,q, namely 


ux = ^P, 


-r = SM?. 


(5.11) 


cty = ^ (cos /3 — cos a) q , 


f3x = -^ (cos a - cos /3) p, 

Py = {sin /3 — sin a) pq , 

qx = ^ (sin a - sin /3) pq . 

The map {X,Y) 1 —)• {t,x) can be constructed as follows. Setting f = x, then / = t in the two 


(5.12) 

(5.13) 


equations at (5.8), we find 

r c = 2 cX^ XX , 

I -c = -2cYx XY , 

respectively. Therefore, using (5.10) we obtain 

XX = 


1 = 2cXxtx, 
1 = -2cYxtY, 


1 

2 X^ 


1 

2 Y^ 


XY = 


1 


_ (l+cosQ:)p 

~ 4 

_ (l+cos/3)ij 

~ 4 

_ (l+cosQ:)p 

~ 4c ’ 


(1+cos^) q 
4c 


(5.14) 


(5.15) 


Given the initial data (2.2), one particular way to assign the corresponding boundary data for 
(5.11)-(5.15) is as follows. In the X-Y plane, consider the line 


70 = {X + T = 0} C 


(5.16) 


parameterized as x >->• (X(x), T(x)) = (x, —x). Along 70 we can assign the boundary data 
{u,a,/3,p,q) by setting 


a = 2 arctan R{0, x), 
j3 = 2 arctan 5(0, x), 


p = l + i22(o,x), 
q = 1 + 52 ( 0 , x). 


u = tto(x), 

at each point (x, —x) G 70 . We recall that, at time t = 0, by (2.2) one has 
ii( 0 ,x) = {ut + c{u)ur,){0,x) = ui{x) + c{uo{x))uo,x{x), 
5(0, x) = (ut - c{u)ux){0,x) = ui{x) - c{uo{x))uo,x{x). 


(5.17) 
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Remark 3. The above construction (5.16)-(5.17) is by no means the unique way to prescribe 
initial values. One should be aware that many distinct solutions to the system (5.11)-(5.15) 
can yield the same solution u = u{t, x) of (2.1)-(2.2). Indeed, let (tt, a, I3,p, q, x, t){X,Y) be one 
particular solution. Let : M i—)• M be two bijections, with cj)' > 0 and "0' > 0. Introduce 
the new independent and dependent variables {X,Y) and {u,a, I3,p,q,x,i) by setting 


X = 0(X), 
{u,aJ,x,i){X,Y) 


Y = V’(T), 
{u,a,P,x,t){X,Y), 


p{X,Y) = p(X,Y)-(^'{X), 


(5.18) 

(5.19) 

(5.20) 


[ q{X,Y) = q{X,Y)-i;'{Y). 

Then, as functions of {X,Y), the variables {u,a, I3,p,q,x,i) provide another solution of the 
same system (5.11)-(5.15). Moreover, by (5.19) the set 


{{i{X,Y),xiX,Y),u{X,Y)y, (X,y)e]R 2 | (5 21) 


coincides with the set in (5.23). Hence it is the graph of the same solution u of ( 2 . 1 ). One 
can regard the variable transformation (5.18) simply as a relabeling of forward and backward 
characteristics, in the solution u. In connection with first order wave equations, relabeling 
symmetries have been studied in [14, 21]. 


Remark 4. The system (5.11)-(5.15) is clearly invariant w.r.t. the addition of an integer 
multiple of 27r to the variables ce, (3. Taking advantage of this property, in the following we 
shall regard a,/3 as points in the quotient manifold T = M/27rZ. As a consequence, we have 
the implications 


a 

/ TT = 

=> COS a 

> 

- 1 

/3 

/ TT = 

=> cos (3 

> 

- 1 


(5.22) 


Remark 5. Since the semilinear system (5.11)-(5.15) has smooth coefficients, for smooth 
initial data all components of the solution remain smooth on the entire X-Y plane. As proved 
in [17], the quadratic terms in (5.13) (containing the product pq) account for transversal wave 
interactions and do not produce hnite time blowup of the variables p, q. Moreover, if the 
values of p, q are uniformly positive and bounded on the line 70 , then they remain uniformly 
positive and bounded on compact sets of the X-Y plane. Throughout this paper, we always 
consider solutions of (5.11)-(5.15) where p,q > 0. 


The main results in [ 8 , 17] can be summarized as 


Theorems. Letc = c{u) be a smooth, uniformly positive function. Let {t,x,u,a, fi,p,q){X,Y) 
be a smooth solution of the semilinear system (5.11)-(5.15) with boundary data as in (5.17). 
Then the function u = u{t, x) whose graph is 

Graph{u) = | (t(A,T), x(A,T), u(A:,T)) ; (A:,y)GM2| (5.23) 

provides the unique conservative solution to the Cauchy problem (2.1)-(2.2). 
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Figure 4: The level sets {a = tt} and {/3 = tt} in a solution with generic singularities. In the X-Y 
plane these are smooth curves which are structurally stable w.r.t. small perturbations. 


Throughout the following, we shall be interested not in a single solution but in a continuous 
path of solutions 6 , 9 £ [0,1], We introduce suitable regularity conditions, allowing 

us to compute the “length” of this path by integrating a suitable norm of its tangent vector 
\\du^{t,-)/de\\. 

Definition 1. We say that a solution u = u{t,x) of (2.1) has generic singnlarities for t G 
[0,T] if it admits a representation of the form (5.23), where (i) the functions {u,a, /3,p,q,x,t){X,Y) 
are C°°, and (ii) on the domain where t{X,Y) G [0,T] the following generic conditions hold: 


(Gl) 

a = TT, 

oix =0 = 

=4> ay / 0, 

axx / 0 

(G2) 

II 

o 

II 

=4- /3x 7^ 0, 

/^VY / 0, 

(G3) 

a = TT, 

Xo 

II 

ii 

» ax 7^ 0, 

/3y 7 ^ 0. 



Figure 5; The set of singular points (where \ux\ +oo) in a solution u{t,x). These are the images 
of the sets {a = tt} and {/3 = tt} in Fig. 4. By structural stability, every small perturbation will yield 
anther solution with the same type of singularities. 

Some words of explanation are in order. Even if the solution {X, Y) i->- (x, t, u, a, f3,p, q){X,Y) 
of the semilinear system (5.11)-(5.15) remains smooth on the entire X-Y plane, the function 
u = u{t,x) in (5.23) can have singularities because the coordinate change A : {X,Y) (x,t) 
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is not smoothly invertible. Indeed, by (5.15)-(5.14), the Jacobian matrix is computed by 


DA 


f XX 

XY \ 

_ 1 

V tx 

ty ) 

■ 1 


(1+cos a) p 
4 

(1+cos a) p 
4c{u) 


(1+cos 0)q 
4 

(1+cos q 
4c{u) 


(5.24) 


We recall that p, q remain uniformly positive and uniformly bounded on compact subsets of 
the X-Y plane. By Remark 3, at a point (Xo,yo) where a ^ tt and (3 ^ tt, this matrix 
is invertible, having a strictly positive determinant. The function u = u{x,t) considered at 
(5.23) is thus smooth on a neighborhood of the point 


{to,xo) = {t{Xo,Yo), x{XQ,Yo)y 

To study the set of points in the x-t plane where u is singular, we thus need to look at points 
where either tc = tt or /3 = vr. The generic conditions (G1)-(G2) guarantee that these level 
sets are smooth curves in the X-Y plane. Gondition (G3) implies that the level sets where 
{a = tt} and {j3 = tt} intersect transversally because ay = I3x = 0 when a = /3 = 0. As 
observed in the conditions (G1)“(G3) are invariant w.r.t. smooth variable transformations 
(A, y) o (A,y). We also remark that, if a solution U = {u, a, I3,p, q) of (5.11)-(5.13) 
satisfies the generic conditions (G1)“(G3), then by the implicit function theorem the same is 
true for every perturbed solution U = {u,a, (3,p,q) sufficiently close to U. In other words, 
generic singularities are structurally stable. An example of structurally unstable solution, 
corresponding to a change of topology in the singular set, is shown in Fig. 6. 


Definition 2. We say that a path of initial data 7 : 0 i—)• {uQ,uf), 9 G [0,1] is a piecewise 
regular path if the following conditions are satisfied. 

(i) There exists a continuous map {X, Y, 6) 1 —)• (tt, a, I3,p, q, x, t) such that, for each 9 G [0,1] 
the semilinear system (5.11)-(5.15) is satisfied. Moreover, the function u^{x,t) whose 
graph is 

Graph{u^) = {(t, x, u)(A, T, 0); (A,y)GM2| (5.25) 

provides the conservative solution of (1.1) with initial data 

u^{0,x) = Uq{x) , u^{0,x) = uf(x). 

fa) There exist finitely many values 0 = 9o < 9i < ■ ■ ■ < 9 n = 1 such that the following 
holds. For 9 G]9i-i,9i[, the map {X,Y,9) (u,a, I3,p,q,x,t) is C°°. Moreover, the 

solution = u^{t,x) has generic singularities at time t = 0. 

In addition, if for all 9 G [0,1] \ {9i,... ,9]\f}, the solution has generic singularities for 
t G [0,T], then we say that the path of solutions ^ \ 9 is piecewise regular for 

t G [o,r]. 


Remark 6. According to Remark 3, there are infinitely many parameterizations of the 
variables (A,T) that yield the same solution u = u{t,x). However, as shown in [7], the 
property of having generic singularities is independent of the particular representation used 
in (5.25). 
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Remark 7. The above definition has a simple motivation. If 7 is a piecewise regular path, then 
we can compute its length as an integral of the norm of a tangent vector. In addition, if 7 is 
piecewise regular fort £ [0, T], then the length of the path of solutions 7 * : 0 1 —)■ {u^{t, ■),uf {t, •)) 
is well defined not only at t = 0 but for all t G [0, T]. See Definition 3 in Section 6 for details. 


Remark 8. In Definition 2, the finitely many values of 6 where does not have structurally 
stable singularities correspond to bifurcation values. As 6 crosses one of these values, the 
topological structure of the singular set (where tt® —>■ ± 00 ) usually changes, as shown in 
Fig. 6 . 





Figure 6 : Here the solution has generic (i.e., structurally stable) singularities ioi 9 < 9 and for 
9 > 9. However, when the parameter 9 crosses the critical value 9, the topology of the singular set 
changes. 

Following [7], on the wave speed c we assume 

(A) The map c : M 1 —)• M+ is smooth and uniformly positive. The quotient c'{u)/c{u) is 
uniformly bounded. Moreover, the following generic condition is satished: 

c'{u) =0 ^ c"{u) / 0. (5.26) 

Notice that, by (5.26), the derivative c'{u) vanishes only at isolated points. The following 
result, proved in [7], shows that the set of piecewise regular paths is dense. 


Theorem 4. Let the wave speed c{u) satisfy the assumptions (A) and letT>0 be given. Let 
6 !-)• {t^,x^,u^,a^, (3^,p^,q^), 9 G [0,1], be a smooth path of solutions to (5.11)-(5.15). Then 
there exists a sequence of paths of solutions 9 1 —)• (t^, x®, a^,< 7 ®) with the following 

properties. 

(i) For each n > 1, the path of corresponding solutions of (2.1) 9 1 —)• is regular for 

t G [0,T], according to Definition 2. 

(a) For any bounded domain D in the X-Y plane, as n —)■ 00 the functions (t®, a®,/3^, 

Pn^Qn) converge to {t^ ,x^ ,u^,a^, (5^,p^,q^) uniformly mC^([0,1] x D), for every k>l. 


Thanks to this density result, to construct a Lipschitz metric it now remains to show that the 
weighted length of a regular path satishes the same estimates as the smooth paths considered 
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in the previous section. Toward this goal, we first derive an expression for the norm of a 
tangent vector as a line integral in the X-Y coordinates. 

Consider a reference solution u (2.1) and a family of perturbed solutions u^, e € [0,eo]- We 
assume that, in the X-Y coordinates, these define a smooth family of solutions of (5.11)-(5.15), 
say ,x^,u^,a^, f3^,p^,q^). For each e, the curves where X =constant and Y = constant 
correspond respectively to backward and forward characteristics of the solutions u^{t,x). We 
remark that, at time t = 0 , we have considerable freedom in choosing these parameterizations. 
We can take advantage of this in the following way. Let w, z be the shifts in (3.26). At time 
t = 0 we choose the parameterizations according to 

A^(0, x + ert;(0,x)) = x, F"^(0, x + ez(0,x)) = —x. (5-27) 

Consider the curve in X-Y space 

r, = {(A,y), t{X,Y) = T} = {(A,y(T,A)); AgM} = {(A(T,y), y); y€M}, 

(5.28) 

and denote by 

= {(A,y), f{X,Y) = T} = {(A,y^(r,A)); AgM} = {(A^(r, A), A); A G M} 

(5.29) 

the perturbed curve. We can write the perturbed solutions as 

{f,x^,u^,a^,/3^,p^,q^) = {t,x,u,a,/3,p,q) s{T, X,U, A, B, P,Q) o{£) (5.30) 

Since the system (5.15)-(5.11) has smooth coefficients, the first order perturbations satisfy a 
linearized system and are well defined for all (A, A) G We observe that the quantities 
v,f,s,w,z appearing in (3.15) can be expressed in terms of the hrst order perturbations 
{T,A;,U,A,B,P,Q). Indeed, 

{1 + R^)dx = pdX , {1 + S'^)dx = — qdY 


Notice that, by dehnition. 


P{X, A^(t,A)) 


t^(A^(r,A),A) 


T. 


Hence by the implicit function theorem, at e = 0: 

dX^ _ dP ( dt\~^ 
de de \9A j 

and 


Qye _ Qte f dt 
de de \9A j 


-T 


4c 


-T 


(1 + cosa)p 
4c 


(1 + cos f5)q 


1. The shift in x is computed by 


w = lim, 


e —>-0 


x'^{X,Y'^{r,X))-x{X,Y{T,X)) 


= A(A,A(r,A)) + xy • 


dYf_ 

de 


£=0 


= (A + cr)(A,A(r,A)). 


In a similar way, 


z = lim. 


£—>0 


x^{X^{t,Y),Y)-x{X{t,Y),Y) 


= A(A(r,A),A) + xx-^ „ = (A - cr)(A(r, A), A), 


£=0 
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2 . We now derive an expression for f, s. One has 


r + wRx = — tan 
de 


d a%X,Y%T,X)) 


= \{A-T 

£=o 2 


4c 


7 -sec — 

(l + cos/3)g / 2 


a 


and 


s + zSx = -r ta-ii 
de 

By (3.14) it thus follows 


d^ f5%X^iT,Y),Y) 


= -(b-t-. -——^ 

£=0 2\ (1 + cosq;)p 


2 P 

sec 2 


and 


1 / . _ 4c 

fp _I /4 _ ' f _ 

2V (l + cos/3)g 


s = -(b-T- - - —- 

2 V (1 + cos a)p 


2 « c' 


ay sec^-T tan^ — 

' 2 4c 2 


2/3 c' 2 « 

sec-— / tan — 


2 4c 


(5.31) 

(5.32) 

(5.33) 

(5.34) 


3. By (5.11) one has 

v + UxW = ■^u^(X,Y^(t,X)) 


^ 4c 

= U — UyT 7 -- 

£=o (l + cosp )(7 


a 


U — T tan — . 


Therefore 


Rw Sz , a /3. ^ 


(5.35) 


4. We now calculate the terms I 4 - Iq in (3.15). 

The change in base measure with density 1 + is given by 
^(X.YHr,Xy,-p(X,Y(r.Xr, ^ 

The change in base measure with density 1 + 5^ is given by 

q%X%T, Y),Y)- q{X{T, Y),Y) ^ , dX^ 

Inn - = Q{X,Y)+qx ■ 

£->■0 e de 


4c 

= P-T --- py . 

£=o (l + cos/3)g 


(5.36) 


4c 


= Q - T- 
£=o (l + cosa)p 


qx 


(5.37) 


The change in base measure with density R? (the integrand in Ig) is estimated by 


^(/sin2^(X,W(r,X))) 
4c 


£ = 0 


= [P-T 


\ 2 ® p sin a / 4c 

(l+COBffl/ d 2 - '^ (l+COBfflP >' ' ■ 


(5.38) 


The difference between (5.36) and (5.38) shows that the change in base measure with density 1 
(the integrand in I 4 ) is computed by 


P-T 


4c \ 9 a p sin a / , ^ 4c 

py cos - - 1—-— A - T 77 —- —ay 


(1 + cos f5)q 


(1 + cos /3)g 


(5.39) 
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Combining the previous computations, the weighted norm of a tangent vector (3.15) can be 
written as a line integral over the line F,- defined at (5.28): 


{f,w,s,z) 


iu,R,S) 


6 

S'"'/ 

e=i 


W- dX + dY 


}> 


(5.40) 


where 

= {X-cT)p 


J2 


h 


Ja 


Jb 




[i+tosP)q ^y) - TcP^ tan2 f cos2 I 
(u — (tan ^ + tan |) -T^p 


Pcos^f 


n 2 

H COS"^ ^ 


psino 

2 


A + ^Tav 


sing 
1+cos 13 


+ tan ^ cos^ ^ 


J 4 •tan j 


P 




' (l+cus 

+ |^(sin^ ^ tan ^ — tan' 


L^fsinf- 


Using (5.13) and (5.12), the above expression can be simplified as 


^ Ji = {X-cT)p 

J2 = ^Ap-^TAn^ 

J 3 = (u — (tan ^ + tan |) • t)p 
< ^ ^ (5.41) 

J 4 = Pcos^ ^ A + I^Tpsina 

J 5 = sin a — pA sin^ f sin^ ^ 

_ Je = Psin^f + P^A. 

In a similar way, we obtain 


Pi 

= (d:' + cr)g 


P 2 



P 3 

= “ (tan § + tan ^ 

)-r)<? 

P 4 

= Qcos^f + 

I^Tgsin/3 

P 5 

= ^Qsin/3 — gPsin^ ^ 

+ fjT q sin- 

Pe 

- gsin^^ + ^^^P 
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It is clear that the integrands Hi are smooth, for I = 1,2,4, 5, 6 . We claim that the 
integrands J 3 and ifa are continuous as well. Indeed, using (5.35) we obtain 

U — (tan ^ + tan ■ T 

= 2cv + Rw — Sz 

= f (^jr(R — S)v + 2cvx + wRx + Rwx — zSx — Szx^ dx 
= f — s + wRx + Rwx — zSx — Szx^ dx 

= f (^r + Rwx - ^{w - 2 ) 5 ^^ dx - f (^s + zSx - ^{w - z)R^'^ dx 
+ f ^wRx + -^R?‘S{w — 2 :)^ dx — J (^Szx + -^R?S{w — z)^ dx 
+ f{£z{w- z){S^ - R‘^)'J dx . 

The three terms on the right hand side correspond to the integrands in I 2 , R and R, respec¬ 
tively. Hence they are continuous. 


6 Length of piecewise regular paths 

Let 7 : 0 I—(uq, tt^) be a piecewise regular path of initial data. According to Definition 2. there 
exists a smooth path of solutions of (5.11)~(5.15), say 6 1 —)• {x^, ,a^, (3^,p^, q^){X, Y), such 

that (5.25) holds for every 9 G [ 0 , 1 ]. At time t = 0 , an upper bound on the length of this 
path can be computed as follows. For each 9 € [0,1], consider the curve in the X-Y plane 

rg = {(x,y); t%x,Y) = oy 

The norm of the tangent vector is then determined by (5.40). Integrating w.r.t. 9 we obtain 
[ [ ^{\Je\^V~dX +\Hf\W+dY^\ d9 . (6.1) 

•^0 \£=1 -^^0 / 

We recall that there exist infinitely many paths of solutions of (5.11)-(5.15) which yields the 
same path of solutions to (2.1). Indeed, as shown in Remark 3, at time t = 0 for each 9 one 
can choose smooth, increasing functions (fP (smoothly depending on 9), and define the 
solutions {x^,F,u^,a^,P^,p^,q^){X,Y) as in (5.18)-(5.20). 

On the other hand, different relabelings of the X, Y coordinates determine different values for 
the integral in (6.1). Indeed, these correspond to different choices of the shifts w,z in (3.13). 
To illustrate this point more clearly, fix a value of 9. Then, for e > 0 small, the family of 
solutions can be regarded as perturbations of the solution . At a given point (r, x), 
the shifts w{t,x) and z{t,x) are uniquely determined as follows (Fig. 7). Let Xq,Yq be the 
point in the X-Y plane such that x^{Xq,Yq) = x, t^{XQ,Y(j) = r. For each e > 0, define 
and Y^ implicitly by setting 

t^^RXo,YR = r, t^^RX„Yo) = T. 
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The shifts are then uniquely defined by setting 

x®+"(Xo,y,)-x^Xo,yo) 


w{t,x,) = lirn 
e^O 


z{t,x) 


lim 
€—^0 


x<^+%X„Yo) - x^{Xo,Yo) 

( 6 . 2 ) 



Figure 7: Given a representation of the solutions in terms of the variables X, T, the shifts w, z are 
uniquely determined by (6.2). Here F® = {{X, Y ); t^{X, Y) = r}. 

The above considerations lead to 


Definition 3. The length || 7 || of the piecewise regular path 7 : 0 i—)• {uq,Ui) is defined as the 
infimum of the expressions in (6.1), taken over all piecewise smooth relabelings of the X-Y 
coordinates. 


Based on the analysis in Section 3, we now give an estimate on how the length of a regular 
path can grow in time. 


Theorem 5. Given any K,T > 0, there exist constants ki, ... ,kq in (6.1) and Ck,t > 0 such 
that the following holds. Consider a path of solutions 6 1 —)• of (1.1), which is piecewise 

regular for t £ [0, T] and where each has total energy < K. Then its length satisfies the 
estimates 


IIt’"!! < C'ty.r IIt'^II for all 0<r<T. 


(6.3) 


Proof. 1. To fix the ideas, let be structurally stable for every 6 £ [0,1] \ {9i ,..., On}- 
Fix e > 0 and choose a relabeling of the variables X, Y such that, at time t = 0, 

+ \Hf\W+ dY^ dO < || 7 °||+e. (6.4) 

Since the solution u is smooth in the X-Y variables and piecewise smooth in the x-t variables, 
the existence of the tangent vector is clear, for every 9 £ [0,1] and t £ [0, T]. We claim that, 
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for every 6 ^ {6\,... , 6 * 7 v}) an estimate such as (3.22) holds. Namely 


(/(r),/(T),s'^(r)) 


{u^{t),R^ (t),S®(t)) 


(6.5) 


< exp<jC'oT+ / a'^(s)ds|>- (n®(0), r'^(O), s^(0)) 


(«®( 0 ),iJ®( 0 ),S®( 0 )) 


Here the constant Cq and the integral of a® depend only on T and on an upper bound on the 
total energy. 

Integrating (6.5) over the interval 9 S [0,1], one obtains an estimate of the form 
h^ll < C(|| 7 °||+e) for all 0<r<T. 

This proves (6.3), because e > 0 was arbitrary. 


2. It now remains to prove the estimate (6.5). We observe that, if were smooth for 
all {x,t) £ M X [0,r], the result follows directly from (3.25), proved by the computations in 
Section 4. We need to show that the same conclusion can be reached if is piecewise smooth, 
with structurally stable singularities. 



Figure 8: Proving that the rate of change in the length of a tangent vector is not affected by the 
presence of a singularity. 

Fix a time r and call F,- = {t^{X,Y) = r} the level set in the X-Y plane. Since the 
estimates of the previous section hold in regions where is smooth, to obtain a bound on the 
weighted norm of the tangent vector it suffices to show that the effect of isolated singularities 
is negligible. To lighten the notation, in the following the superscript ® will be omitted. 

With reference to Fig. 8 , assume that the solution has a structurally stable singularity along 
a backward characteristic. We claim that this singularity does not affect the estimate (3.25). 
In other words, the time derivative 

,6 

/ {\M^~ dX + dY^ 

is not affected by the presence of the singularity. 
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Figure 9: Here P 2 is a singularity point of Type 2, where a = tt and ax = 0, but axx 7 ^ 0 and 
j3 ^ TT. At P 3 the solution has a singularity of Type 3, where a = /3 = tt, but ax 7 ^ 0 and Py 7 ^ 0. The 
weighted norm of the tangent vector is continuous at the times t = t 2 and t = 


For a given time r, let be the point where the curve Tr-e = = t — e] 

intersects the singular curve {a{X,Y) = tt }. Similarly, let (X'^,Y^) be the point where the 
curve Ft+s = {t(A", Y) = t + e} intersects the singular curve {a{X, Y) = vr}. 

Define the curves 

r C7+ = F,+, n {X e [X',X,]}, r r/+ A n {y g [y„y,]}, 

1 u- A n {x G [X', X,]}, I r/r A n {y g [y„ y,]}. 

To prove our claim, it suffices to show that 


lim — ( 

£— >-0 S \ 

-..JJUq J 

( 6 . 6 ) 

lim — f 

£—>-0 S \ 

f - f ) = 0- 

Pvt drje J 

(6.7) 


The first limit holds because the integrand is a continuous function of X, y and \X^ — 
X'l = 0(e). The second limit holds because the integrand is a continuous function of 

X, y and \Y^ — Ye\ = 0(e). The basic estimate (3.25) thus remains valid also in the presence 
of singular curves where a = n ox (3 = tt. 


Finally, we analyze what happens in the presence of singular points of Type 2, where a = vr 
and ax = 0, and of Type 3, where a = P = tt. Since the solution is structurally stable, 
there can be at most finitely many such points, say 

Q, = (x,-,y,), j = i,...,x. 

To complete the proof of our claim, it thus suffices to show that, at each time Tj = t(Xj,Yj), 
the map 

^ ^ J {\Ji\^~dX + \Hf\W+dY^'^ ( 6 . 8 ) 
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is continuous at t = Tj. But this is clear, because the path depends continuously on t and 
the integrands Ji, Hi are uniformly bounded. Moreover, they are continuous everywhere with 
a possible exception of the finitely many singular points Qj. □ 


7 Construction of the geodesic distance 


A key result proved in [7] shows that every path of solutions to (1.1) can be approximated 
by a path which remains regular for t G [0,7"]. More precisely, an application of Thom’s 
transversality theorem yields 


Theorem 6. Let the wave speed c{u) satisfy the assumptions (A). Let (u®, ,p^,q^,x^, Y) 

be a path of C°^ solutions to the semilinear system (5.11)-(5.15), depending smoothly on 
6 G [0,1]. Then, for any T,s > 0 and any integer k > 1, there exists a perturbed path of 
solutions {vf ,oP,pP,x^,?){X,Y) such that 


{u^ - - a\ / 3 ^ - / 3 ^ / - /, q^ - q\ x^ - x\ t^ - ?) 

Here C is a domain containing the set 


c^{n) 


< e. 


(7.1) 


|(A:,T); t®(X,y) G [0,r] or i\X,Y) /or some 0 G [0,1]I. 

Moreover, all except finitely many solutions {vP,p^,x^,i^) have structurally stable 
singularities inside tl. 


In other words, by slightly perturbing the initial data (uq,u^), 6 G [0,1], we can construct 
a one-parameter family of conservative solutions = u^{t,x) which have structurally stable 
singularities, for all but finitely many values of 6 . This implies that for all t G [0,T] the length 
of the path 6 (t, ■) is well defined by the formula 


de 


«»(() 


M. 


U^(t) 


(7.2) 


Here || • ||„ is a weighted norm defined as in (3.13)-(3.15), or equivalently at (5.40). 


A geodesic distance d* on the space H^(W) x L^(M) will be constructed in two steps. 

(i) As proved in [7], there is an open dense set of initial data 

V c (c^{R)nH^{R)'^ X (c^{R)nL‘^{R)y (7.3) 

such that, if {uo,ui) G V, then the solution of (2.1)-(2.2) has structurally stable singu¬ 
larities. On n P we construct a geodesic distance, defined as the infimum 

among the weighted lengths of all piecewise regular paths connecting two given points. 

(ii) By continuity, this distance can then be extended from to a larger space, defined as 

the completion of w.r.t. the distance d*. In particular, this completion will contain 
the space {H^ 0 x (L^ n L^). 
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More in detail, assume {uq,ui), {uq,ui) G V°°. Their total energies will be denoted by 

£{uo,ui) = J [ul + c^{uo)ul ,,]dx , £{uo,ui) = J [ul + c^{uo)ul ,,] dx, 

respectively. Fix any constant K > 0 and consider the subset of all data with energy < K, 
namely 

Xk = {{uo,ui)gH\R)xL\R)-, £{uo,m)<KY (7.4) 

Notice that Xk is positively invariant for the flow generated by the wave equation. 


Definition 4. On 7)°° Ci Xk we define the geodesic distance d*((uo,ui), {uQ,ui)'j as the 
infimum among all weighted lengths of piecewise regular paths, which connect {uo,ui) with 
(no,iii), always remaining inside Xk- Namely, 


d*((uo,ui), (uo,ui)') = inf |||7||; 'y is a piecewise regular path, 

7(0) = (rio,rii), 7 ( 1 ) = (fto, %) £(ul,u\) < K for all 6 e . 


(7.5) 


Since the concatenation of two piecewise regular paths is still a piecewise regular path (after 
a suitable re-parameterization), it is clear that d*{-,-) is indeed a distance. As a consequence 
of Theorem 5, we have 


Theorem 7. Let the wave speed c(-) be smooth and satisfy (2.3). Then the geodesic distance 
d* renders Lipschitz continuous the flow generated by the wave equation (2.1). In particular, 
let {uo,ui) and {uo,ui) be two initial data in (2.2). Then for all t G [0,T] the corresponding 
solutions satisfy 

d*(^{u{t,-), uflt,-)), {u{t,-), ut{t,-))^ < Ck,t ■ d*(j^uo,ui), {uo,ui)J . (7.6) 

Here Ck,t is a constant depending only on T and on an upper bound K on the total energy. 

Proof. If the wave speed c(-) satisfies the generic assumption (A) at (5.26), then the result is 
a direct consequence of Theorem 5. To cover the general case, it suffices to approximate c(-) 
with a sequence of functions c„(-) that satisfy the assumption (A). If \\cn — c||c3(q) —)• 0 as 
n —)• 00 for every bounded interval C M, then the flow generated by the velocities c„(-) and 
the corresponding geodesic distances converge to the ones for c(-). □ 


In the remainder of this section we compare the distance d* with more familiar distances in 
Sobolev spaces, and with a Wasserstein distance between energy measures. 

Proposition 2. There exists a constant C'^ such that, for any (uo,ui), {uo,ui) G 'D°°r\XK, 
d*[{uo,ui), {uo,ui)) < (||^to-'Uo||Hl + ll^^o-'Wo||M/l■l + ll^^l-'“l||L2 + ||^^l-'Ml||Ll)• (7.7) 
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Proof. 1. Define the function 


pu 

'^{u) = / c{s) ds. (7.8) 

Jo 

Observe that : M i—>■ M is a smooth strictly increasing function, with smooth inverse '1'“^. 
The total energy can then be expressed as 


£{uo,ui) = J[uj + c^{uo)ul^^]dx = j[u‘l +{^{uo)x)‘^]dx. 


Let (ho, hi) be another initial data, with total energy £. For 6 G [0,1], consider the interpolated 
data {uQ,uf) where 


= T 


-1 


(0T(ho) + (l-0)T(uo)), 


(7.9) 


9ui + (1 — 9)ui. 


When 6 = 0,1, it is clear that (ttQ,ttf) coincides with (rto,tti) and (ho, hi), respectively. We 
check that the energy remains < M. Indeed, 


j [{Uif + c^(«o)K,x)^] + (^K)^)^] dx 

= J [{6ui + {1 - e)uif dx + I [9^{uo)x + {I - 9)'^{uo)x]^ dx 
< max {T(uo, ui), T(ho, hi)} < M. 


2. Next, we estimate the weighted length of the path 7 : 0 i—)• {uQ,uf) in (7.9), showing that 


< C ■ M|no - holl^i + Iko - 'Wollwm + ||'«i - ui\\i^2 + 


(7.11) 


for some constant C depending only on the total energy. To establish an upper bound for the 
weighted length || 7 ||, in the definition (3.15) we choose the shifts w = z = 0. In this way, the 
integrals and I 5 vanish. 


We first calculate {v^,r^,s^) = ,R^,S^). From (7.8) it follows 

'^'{u) = c{u), 

(«)) ^'(^-i(a)) c(T-i(a)) ' 

Using (7.9) and (7.12)-(7.13) we find 


^ T(ho) - T(uo) 

d9 c(0'I'(ho) + (1 - 0)'I'(tto)) 


Since the wave speed c(-) is uniformly positive, the above implies 


Ki 


|ho-no| < |/| < iLi|ho-tto|, 


(7.12) 

(7.13) 


(7.14) 


for a suitable constant Ki, depending on the function c(-) and on an upper bound for the 
energy. 
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Next, we have 


= u\ + ^{4), = 0(ui + ^(no)x)+ (l-0)(ui + ^'(uo)x) = eR+il-e)R. (7.15) 

Hence 


d 


= —R^ = (mi + ^('Uo)a;) - (ui + ^(ito)x) = R-R. 


Similarly, 


= {ui - 'h(no)x) - - 'h(^xo)x) = s-s. 


(7.16) 

(7.17) 


For later use, we observe that 


^ (^j\2R^r^\dx^d9 = j\R-R\-(^j^ 2\R^\de^dx< j\R-R\ ■ {\R\+ \R\) dx . (7.18) 

Observing that the weights satisfy a uniform bound depending only on the total energy, 
and using (7.14)-(7.18), we finally obtain 


{uO,Rf>,S^) 


r {v\r\s^) 

Jo 

J [\r^\{W-f + \s^\{W+f]dx 

+K3 J\v^\{{l + (Rff) {W-f + (1 + (Sy) (W+)®} dx 
+kq j ||2i?V|(W-)® + |25^s®|(W+)'^}dxld0 


< K2 


.R - ii| -I- |5 - S'! \dx + ||uo - uoIIli 


-l-||no — fio||L°° 


^ ( j {(R'^)^ + (R®)^}(ix)d0 


(7.19) 


+ / |R - R| • (|R| + |R|) + |S - 5| • (|5| + |S|) dx 


< ■ (||uo — UqWhi + ll^^o — 'dolIvKi'i + 11^1 “ 'dl||L2 + 11^1 — 'WiIIli)) 

where K 2 and K 3 are positive constants, depending on the upper bound for the energy. In 
the last step, we use similar estimates as in (7.10). This completes the proof. □ 

We conclude the paper by showing that the geodesic distance d* in (1.4) controls both the 
distance ||uo —uqIIli the Wasserstein distance between the corresponding energy measures 

/r,/l. 
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Proposition 3. There exists a constant 6 q, depending only on an upper bound on the energy, 
such that for any uq, uq S n and any ui,ui € one has 


|| wo - wo||li < dQ-d*{{uo,ui), {uQ,ui)), 


(7.20) 


sup 

ll/llci<l 



< (io • d*{{uQ,ui), {uQ,ui)) . 


(7.21) 


Here pL,ft are the measures with densities u\ + c^(uo)uoa: + c^(^o)'^oa: w.r.t. Lebesgue 

measure. 


Proof. 1. To prove (7.20) we first observe that 


|u| < 


V + 


Rw Sz 
2 c 2c 


+ 


Rw 


2 c 


+ 


Sz 


2 c 


< 


V + 


Rw Sz 
2 c 2c 




(7.22) 

The right hand side of (7.22) is bounded by the integrands in R and R in (3.15). Recalling 
the definition (7.5), by (7.14) for some constant C 4 > 0 we thus have 


d*{{uo,ui), {uo,ui)) > C4-inf 



\v^ \ dx dO 


< C 4 • inf 

7 


f 


du^ 


(7.23) 


de 


d9 = C 4 ||uo — uqIIli 


2 . Next, consider any regular path 7 : 0 1 —)• (ttg,u®) joining (uo,ui) with {uo,ui). Call pP the 
measure having density {uf)^ + c^(nQ)(uQ)^ = (R®)^ + (R^)^ w.r.t. Lebesgue measure. 

Then, for any function / such that ||/||ci < Ij one has 



< f \f\-{H{^ + R^) + \^\{^ + s^)}dx 


+R5 

• ^ I/I • 2 R{r + wRx) + Rpwx + 2 S{s + zSx) + S'^Zx) 


< Rs- 

y 1 w (1 + R^) + \z\{l + S^)'^dx 


+R5 

■ j 1 2 R(r + wRx) + R^Wx + “ S‘^R){w - z) 



(7.24) 



dx. 


Using (3.14), we see that the two integrals on the right hand side of (7.24) are exactly R 
and Iq without potential terms LV~ and W'*', hence are dominated by the integrals in (3.15). 
Integrating w.r.t. 9 G [0,1], one obtains (7.21). □ 
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